Large Angle CMB Fluctuations from Cosmic Strings with a Comological Constant 



cn 

O 

o 



(N 
> 

(N 

O . 
m . 
o ■ 

i!^ '■ 

Oh 
I 

o 



M. LandriaiQ 

Laboratoire de I'Accelerateur Lineaire, 
IN2P3-CNRS et Umversite Pans-Sud, 
B.P. 34, 91898 Orsay Cedex, France 

E.P.S. Shellarcill 

Department of Applied Mathematics and Theoretical Physics, 
Centre for Mathematical Sciences, 
University of Cambridge 
Wilberforce Road, Cambridge CB3 OWA, U.K. 
(Dated: August 1, 2003) 

In this paper, we present results for large-angle CMB anisotropies generated from high resolution 
simulations of cosmic string networks in a range of flat FRW universes with a cosmological constant. 
Using an ensemble of all-sky maps, we compare with the COBE data to infer a normalization (or 
upper bound) on the string linear energy density p. For a flat matter-dominated model {Q,m ~ 1) 
we find Gfi/c^ ~ 0.7 x I0~^, which is lower than previous constraints probably because of the 
more accurate inclusion of string small-scale structure. For a cosmological constant within an 
observationally acceptable range, we find a relatively weak dependence with G^,/c? less than 10% 
higher. 

PACS numbers: 98.80.-k, 98.80.Cq 



I. INTRODUCTION 

Given the growing evidence for a cosmological constant 
or dark energy component in the Universe, it is impor- 
tant to determine whether this significantly impacts the 
evolution of cosmic strings and their observational signa- 
tures. The strongest constraints limiting the energy scale 
of cosmic strings is currently the COBE normalisation of 
large-angle CMB anisotropy. In this paper, we will infer 
this normalisation from maps computed using the meth- 
ods presented in I'l and determine how this is influenced 
by late-time domination by a cosmological constant. 



II. METHOD 

A. Geometry of Simulations 

To compute all-sky maps, the "observers" are located 
inside the simulation box. Following 0, we place eight 
such observers at each corner of a cube of side L/2 where 
L is the size of the box itself. We use the fact that the 
string networks, and hence the cosmological perturba- 
tions they induce, have periodic boundary conditions to 
identify opposite sides of the box. The eight realizations 
of the sky are not completely independent, but the corre- 
lations are not apparent. This scheme reduces the effect 
of the cosmic variance on large angular scales and thus 



renders the computation of the normalization more ac- 
curate. 



B. Pixelization 

To facilitate the computation of maps and the extrac- 
tion of the power spectrum, we use an iso-latitude pix- 
elization, which enables us to Fourier transform each ring 
of constant 9 to evaluate the integral over the azimuthal 
angle. More specifically, we use the equal area with 3:6:3 
base pixels, proposed by 0]- In this scheme, each base 
pixel is subdivided into = 2"^^ pixels. The area of a 
pixel is given by Api^ei = ■^'^/Np^xei = Tr/Sn^. But, by 
integrating the spherical surface element over a pixel, we 
have Apixei = -^AcosS. Hence, the width of a row of 
pixels is given by 



Acos0 = ^, 



(1) 
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where Np is the number of pixels per row of constant 9. 
In the polar caps, Np — 3, then 9 x 2*~^ repeated 2*~^ 
times with 2 < i < p + 2. The equatorial regions are 
made up of n rows of 6n pixels. 



C. Angular Power Spectrum Computations 

The pixelized map is the convolution of the real tem- 
perature map with the pixel window function. To com- 
pute the power spectrum, we need to know the spherical 
transform of the window function. In real space, the pixel 
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window function is simply given by: 
' 1 outside 



(2) 



where p and q label the rows and the pixels on that row: 

(3) 



TP? = / T{<p,e)WP'i{(j),9)dn 

IS2 



where 



f=|m 



with 



■t. COS f/p — COS Wp+l J COS Op ^ ^ ^ ^ 



(4) 



(5) 



From the last line in J^Jl, the hm can be obtained by the 
inverse sum: 



-imq2'K /Np 



(6) 



which can be done using a Fast Fourier Transform. The 
sum (0J) can also be inverted to obtain the a^m's: 



X;p=T° w*{TTm/Np)w{Trm/Np)Np 



(7) 



To obtain this last equation, we took the approximation 
that the pixelized spherical harmonics w{'Krn/Np)X'^ are 
orthogonal, which is valid at low £. 



D. COBE Normalisation 

To infer the cosmic string linear energy density, we use 
the COBE angular correlation function because, on the 
scales of interests, it is good as that of WMAP This 
also enables direct comparisons with previous work. 

To normalize to COBE, we follow 5] and compute the 
angular correlation function using 

C{9,9s) = ^;^(2^+1)|M^,|2|B,|2qP,(cos0), (8) 



where we use the following smoothing function: 



We{0) = exp 



log 2 



sin2(6l/4) 



(9) 



where 6 is the smoothing scales. We also use the COBE 
beam profile Bi given in ^Qj.We then compare the result 
with the value obtained by the four year COBE maps 0, 
C(0°, 10°) = (29 ± IfiK/Tf to infer the value of G^i/c^. 
In the previous equation, T is the mean temperature of 
the CMB also obtained from COBE H T = 2.726 ± 
0.002i^. 



III. RESULTS 
A. Cosmic String Simulations 

We created a time series of cosmic string realizations 
using the AUen-Shellard network simulation 9] in a vari- 
ety of flat FRW universes with il„i + f^A = 1, where flm 
includes the contribution of dark matter and baryons. 
These string networks were then used as sources for the 
CMB map-making pipeline. For all the simulations, we 
chose h = 0.72 and Qth^ = 0.02. 

Two string simulations for JIa = and JIa = 0.7 were 
computed with very high resolution with over 3 million 
points (16 ponts per correlation length) and spanning a 
dynamic range of 5 in conformal time, probing the Uni- 
verse back to z = 25 and z = 30 respectively. All string 
loops were retained in the simulation to ensure overall 
network energy-momentum conservation, which was pre- 
served at the level of 1%. Point-joining was used only to 
ensure a consistent physical resolution, i.e. removing very 
short segments formed through reconnection. The simu- 
lations therefore retained their full small-scale structure. 
These were the costliest part of the pipeline and took 
many days to perform on the COSMOS supercomputer 
(approximately 2000 hours of CPU time). This improves 
upon previous work with string networks in an expanding 
Universe, which was performed with fewer points per cor- 
relation length, fixed horizon resolution and a separate 
treatment of loops 0, ^ . 

In addition, we performed a series of smaller 
simulations (1 million points) for different values 
of the cosmological constants, specifically, JIa = 
0, 0.2, 0.4, 0.55, 0.7, 0.85. These simulations had iden- 
tical initial conditions, so comparisons could minimize 
the effects of cosmic variance. 

The map-making pipeline was tested at resolutions 
of 128'^ and 256"^. The huge data storage require- 
ments with string networks, SVT-decomposed energy- 
momentum tensor grids, Boltzmann code output plus 
checkpointing files meant disk space approaching 1Tbyte 
was required for 256^^. Further parallelisation issues re- 
main to be resolved for large memory runs. 

In figure ^ we present realizations for the fl\ — 
and 0.7 runs. The angular power spectra for each of 
these runs are presented in figure 01 We also show four 
spectra obtained from the ^Ia — run with their average 
in figure [5] All these spectra are normalised to COBE. 



B. Normalization to COBE 

The power spectra illustrated in figures |21 and |31 show a 
roughly scale invariant plateau for £ <20 with a gradual 
fall-off to an insignificant signal beyond £ > 50. This 
fall-off is a consequence of the limited dynamic range of 
the string simulations and not of the resolution at which 
the Boltzmann evolution was performed (though there is 
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FIG. 1: All-sky string-induced maps for a flat FRW Universe 
with f^A = and. 0.7 




FIG. 2: Individual power spectra for the JIa — run in blue 
with the average in red. 



some influence from the interpolation and smoothing of 
the string network onto the grids. 

Previous CMB work with string networks over a much 
larger dynamic range has demonstrated that the main 
anisotropics relevant to COBE are generated at redshifts 
z < 20 0, which is within the range of the present 
work. The reason for this can be seen from studies of 
the unequal-time correlators (UETC's) [H [l3|, which 
peak on scales well below the horizon {Xmax < dn/S). 
We are confident therefore that our simulations include 
the primary contributions on COBE scales. 

This scale invariance on large angular scales is con- 
sistent with the COBE data and allows us to normalize 
our results accordingly, as discussed above, providing a 
constraint on the energy density of strings. Our normal- 





FIG. 3: Average power spectra for each run. 



ization for the flat CDM model. 



^ = (0.7 ±0.2) X 10- 



(10) 



has an error reflecting the variance of our results, not 
systematic effects which may be comparable. This result 
is lower but consistent with the previous result |i5| : 



% = (i.05t°:f ) X 10 



(11) 



Our result is significantly lower than the one obtained 
in Gn/c^ = 1.7 X 10~^ from a single simulation. 
This work however focused on small angle anisotropics 
and did not study sample variance at large angles. Other 
fiat space approximations and semi-analytic estimates for 
local string networks generally have found a higher nor- 
malization than ours, e.g. Gul/c^ = (1.5±0.5) X 10~^ [l3|. 
G/i/c2 = (1.7 ± 0.7) X 10-6 [g and G^/c^ = 2 [13. 

However, there is a key reason why the present work is 
a significant advance over these previous analyses. Apart 
from now including all the relevant physics, this analy- 
sis uses the highest resolution string simulations to date. 
The initial points per correlation length have been cho- 
sen at levels known to preserve small scale structure ac- 
curately ^3 and the simulations are evolved at fixed 
physical resolution. Although the simulations used in 
spanned a larger dynamic range, z < 100 (compared 
to our z < 25), it did so by maintaining fixed horizon 
resolution through point-joining and smoothing. Sub- 
sequent work has shown that at least comoving resolu- 
tion is required if we are to hope for a satisfactory treat- 
ment of string wiggiiness. This implies that the previous 
work did not adequately account for the renormalized 
string energy per unit length, which in the matter era 
is /i ~ 1.4/z. Such a factor would tend to increase the 
string anisotropy, thus lowering the COBE normalization 
of G/i/c^ (though in a non-trivial manner). 
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FIG. 4: Average rms long string velocity for the last half of 
the string network simulations for Qa ~ (red), nA = 0.4 
(green), SIa = 0.7 (blue) and Ha = 0.85 (black). Even the 
extremal A-model shows only a relatively small 11% decline 
in the average string velocity by the present day. (Relative 
conformal time is plotted with rjo = 0.5.) 

C. Effect of the Cosmological Constant 

The influence of JIa on the COBE normalization, illus- 
trated in figureElis relatively small: for the popular value 
of = 0.7 we obtained G^i/c^ = (0.74 ± 0.20) x 10-^ 
only 6% higher than for the flat CDM model. 

The reason for the reduced anisotropy in A-models is 
fairly clear. As the Universe becomes vacuum dominated 
at late times, the expansion rate increases, affecting the 
Hubble damping term a/ a in the string equations of mo- 
tion, thus lowering their average velocity as can be seen 
from figure 21 As string velocities are the primary cause 
of CMB anisotropies, there will be a net reduction in 
AT/T. 

The smallness of this effect is somewhat surprising, but 
is explained by the late redshift of vacuum domination. 
In the r^A — 0.7, this occurs at za — 0.33, which, in the 
context of our simulations, implies that A has a signif- 
icant effect on only the later stages of the simulation. 
Comparing the SIa — 0.85 model with the flat CDM 
model in figure O there appears to be a slight relative 
fall-off in the average power towards i = 2, which is con- 
sistent with this picture. However, this effect is likely to 
be swamped by cosmic variance. 

We have found that our data points are well fitted by 

^i. („..., ^)..„-« 

This analytic fit is motivated by the asymptotic limit in 
which string velocities should vanish as the network is 
frozen and conformally streched in an extreme A-model: 
r^A — > 1 and flfn — > 0. Our results are summarized in 
tableland plotted in figure |S1 along with the fit H12I) . 

IV. SUMMARY 

In this paper, we presented large angle maps of CMB 
fluctuations seeded by networks of cosmic strings in flat 
FRW universes with a cosmological constant. From the 



Gfi/c^ X 10"'^ 

0.00 0.7038 ± 0.1947 

0.20 0.7170 ± 0.1926 

0.40 0.7133 ± 0.1989 

0.55 0.7189 ± 0.1902 

0.70 0.7389 ± 0.1990 

0.85 0.7766 ± 0.1879 



TABLE I: String linear energy density obtained from the 
COBE normalisation. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 5: String linear energy density obtained from the COBE 
normalisation. The solid line is the curve l|12|l. 



COBE data, we have obtained a constraint on the string 
linear energy density G/i/c^ < 0.7 x 10~® which is lower 
than previous work because string small-scale structure is 
incorporated more accurately in the network simulations. 
We were able to find the dependence of the string density 
/i as a function of the cosmological constant, obtaining a 
good fit with a simple semi-analytic formula. Given the 
uncertainties in the overall normalization, these results 
should provide an adequate means by which to charac- 
terise the effects of flA on cosmic string models. 

Although current CMB data have shown that cos- 
mic strings cannot be the dominating source of CMB 
anisotropies (e.g. they can nonetheless be present 

albeit at a lower energy scale. For example, they are co- 
piously produced at the end of brane inflation 0, . 
With this in mind, our normalisation appears low when 
compared with the value inferred from the possible detec- 
tion of a gravitational lensing event by a cosmic string: 
Gfi/c^ w 0.4 X 10~^ 113. This hnear mass density would 
be somewhat higher than allowed in the standard cosmic 
string scenario considered in this paper. 



5 



Acknowledgements 

We are grateful for useful discussions with Gareth 
Amery, Richard Battye, Martin Bucher, Carlos Martins, 
Proty Wu and Rob Crittenden. The AUen-Shellard string 



simulation was used to generate the networks used as 
sources in this paper fs|. This work was supported by 
PPARC grant no. PPA/G/O/1999/00603. All simula- 
tions were performed on COSMOS, the Origin 3800 su- 
percomputer, funded by SGI, HEFCE and PPARC. 



[1] M. Landriau and E. P. S. Shellard, Physical Review D67, 
103512 (2003). 

[2] U. L. Pen, D. N. Spergel, and N. Turok, Phys.Rev. D49, 
692 (1994). 

[3] R. G. Crittenden and N. G. Turok (1998), astro- 

ph/9806374. 
[4] M. Halpern (2003), private communication. 
[5] B. Allen, R. R. Caldwell, E. P. S. Shellard, A. Stebbins, 

and S. Veeraraghavan, Phys.Rev.Lett. 77, 3061 (1996). 
[6] E. L. Wright et al., ApJ 420, 1 (1994). 
[7] A. J. Banday et al., ApJ 475, 393 (1997). 
[8] C. L. Bennett et al., ApJ Letters 464, LI (1996). 
[9] B. Allen and E. P. S. Shellard, Phys.Rev.Lett. 64, 119 

(1990). 

[10] B. Allen, R. R. Caldwell, S. Dodelson, L. Knox, E. P. S. 
Shellard, and A. Stebbins, Phys.Rev.Lett. 79, 2624 
(1997). 

[11] A. Albrecht, D. Coulson, P. G. Ferreira, and J. Magueijo, 



Phys.Rev.Lett. 76, 1413 (1996). 

[12] J. H. P. Wu, P. P. Avelino, E. P. S. Shellard, and B. Allen, 
Int.J.Mod.Phys. Dll, 61 (2002). 

[13] D. P. Bennett, A. Stebbins, and F. R. Bouchet, ApJ Let- 
ters 399, L5 (1992). 

[14] L. Perivolaropoulos, Phys.Lett. B298, 305 (1993). 

[15] D. Coulson, P. Ferreira, P. Graham, and N. Turok, Na- 
ture 368, 27 (1994). 

[16] C. J. A. P. Martins and E. P. S. Shellard (2001), unpub- 
lished. 

[17] C. L. Bennett et al., ApJS 148, 1 (2003). 
[18] S. Sarangi and S. H. H. Tye, Phys. Lett. B536, 185 
(2002). 

[19] L. Pogosian, S. H. H. Tye, I. Wasserman, and M. Wyman, 

Phys. Rev. D68, 023506 (2003). 
[20] M. Sazhin et al. (2003), astro-ph/0302547. 



